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An earlier calculation in a generalized linear sigma model showed that the well-known current 
algebra formula for low energy pion pion scattering held even though the massless Nambu Goldstone 
pion contained a small admixture of a two-quark two-antiquark field. Here we turn on the pion mass 
and note that the current algebra formula no longer holds exactly. We discuss this small deviation 
- - . and also study the effects of an SU(3) symmetric quark mass type term on the masses and mixings 

of the eight SU(3) multiplets in the model. We calculate the s wave scattering lengths, including the 
I beyond current algebra theorem corrections due to the scalar mesons, and observe that the model 

>— ^ , can fit the data well. In the process, we uncover the way in which linear sigma models give controlled 

' corrections (due to the presence of scalar mesons) to the current algebra scattering formula. Such 

{^JQ' a feature is commonly thought to exist only in the non-linear sigma model approach. 

^ : 

^ . PACS numbers: 13.75.Lb, ll.15.Pg, ll.80.Et, 12.39.Fe 

^ : 

. I. INTRODUCTION 

' A linear sigma model with both quark-antiquark type fields and fields containing (in an unspecified configuration) 
two quarks and two antiquarks, seems useful for understanding the light scalar spectrum of QCD. In a previous 
treatment [l| we considered a usual simplification in which the three light quark masses were taken to be zero. The 
1^ ■ model was seen to give a neat intuitive explanation of how "four quark" scalar states could be naturally much lighter 
i-S^ than the conventional p-wave quark-antiquark scalars. We also verified in detail that, as long as the potential of the 
model satisfied SU(3)l x SU(3)r invariance, the massless version of the famous current algebra theorem [2] on low 
energy pion pion scattering was correct. 
^ ' In the present paper we introduce a common mass for the three light quarks in such a way that the pion gets its 
, correct mass. SU(3) fiavor invariance continues to hold, which is a desirable simplification. First we reexamine the 
CD 1 masses and mixings of the particles in the model. It is seen that the natural explanation for the lightness of a "four 
quark" scalar remains unchanged. Then we reexamine the pion pion scattering amplitude to try to see if the low 
energy theorem continues to hold. Curiously we find that it does not exactly hold. What goes wrong? The algebra of 
OO ' the Noether currents should be good in this model so that is not the cause. It turns out that the partially conserved 
axial vector current, which is also required for the theorem, does not hold, unlike for the massless case. The axial 
vector current has a single particle contribution from the "heavy pion" in this chiral model as well as from the ordinary 
pion. The ordinary pion does not therefore completely saturate the axial current. Actually this is a small effect but 
is of conceptual interest and may be of more importance for the kaon scattering case. 

A more important quantitative effect arises from the contributions of the scalar isosinglet mesons in the model. It 
is shown that these contributions can explain the experimental s wave isosinglet scattering length. 
d The notation is reviewed in Section II. The corrections to the masses and mixings, due to the quark mass term, 

are studied in Section III. An approximate analytic treatment of the scattering, for a general potential, is contained 
in Section IV. The exact numerical treatment for a "leading order" potential is presented in Section V. Some discus- 
sion and conclusions are given in Section VI. The Appendix explains the method of parameter determination from 
experiment and a listing of typical values for all the parameters. 
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II. NOTATION 



We introduce the 3x3 matrix chiral nonet fields; 

M^S + icI), M'^S' + i(l)'. (1) 

Here M represents scalar, S and pseudoscalar, quark-antiquark type states, while M' represents states which are 
made of two quarks and two antiquarks. The transformation properties under SU(3)lX SU(3)rX U(1)a are 

M e^"" UlMUI, M' e-'^"' UlM'U^, (2) 

where Uh and C/r are unitary unimodular matrices, and the phase v is associated with the U(1)a transformation. 
The general Lagrangian density which defines our model is 

C = -^Tr (d^Md^M^) - ^Tr {d,,M'd^A'I'^) - Vo {M, M') - Vsb, (3) 

where Vo{M, M') stands for a function made from SU(3)lX SU(3)r (but not necessarily U(1)a) invariants formed out 
of M and M' . The quantity Vsb stands for chiral symmetry breaking terms which transform in the same way as the 
quark mass terms in the fundamental QCD Lagrangian. In our previous paper we focused on general properties 
which continued to hold when Vsb was set to zero. Here, we include the SU(3) symmetric mass term: 

VsB^~2ATr{S) (4) 

where A is a real parameter. A characteristic feature of the model is the presence of "two-quark" and "four-quark" 
condensates: 

{Sl)=aa5l {S',')^f3aSl (5) 

We shall assume the vacuum to be SU(3)v invariant, which implies 

ai = a2 = as = a, Pi = f32 = = P- (6) 

The SU(3) particle content of the model consists of two pseudoscalar octets, two pseudoscalar singlets, two scalar octets 
and two scalar singlets. This gives us eight different masses and four mixing angles. We next give the notations for 
resolving the nonets into SU(3) octets and singlets. Note the matrix convention (ffi^ — + (f)ab- The properly normalized 
singlet states are: 

00 = ^Tr(0), cp'^ = ^Tr(0'), 



Then we have the matrix decompositions: 



7!'^"^' '^' = '^' + 7! 



S^S + ^SqI, S'^S' + ^S'^l, (8) 
V3 \/3 

wherein 0, 0', S and S' are all 3 x 3 traceless matrices. The singlet scalar fields may be further decomposed as: 

5o - VSa + So, S'q^V?>I3 + S'o. (9) 

Here So and Sq are the fluctuation fields around the true ground state of the model. The breaking of SU(3) to the 
isospin group SU(2) will be examined in the future. In that case there are 16 different masses, four 2x2 mixing 
matrices and two 4x4 mixing matrices. To fully characterize the system we will also require some knowledge of the 
axial vector and vector currents obtained by Noether's method: 

(Jr'na = i{c^a-ab)d^S'^ + i{Pa-f3t)d^S'l + ---, (10) 
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where the dots stand for terms biHnear in the fields. 

In our model we use a previously discussed scheme to select the most important terms in the potential, Vo{M, M') 
The favored terms which are SU(3)lxSU(3)r invariant but violate U(1)a are: 



= C3 [F^{M,M')]\ 



in which C3 is a coupling constant and 



Fr,iM,M')=^i In 



/ det(Af) \ 
Uet(Aft) J 



(l-7i)ln 



/Tr(MM't) 
lyTr(M'Mt) 



(11) 



(12) 



where 71 is a dimensionless parameter. This form exactly mocks up the U(1)a anomaly of QCD. Information about 
the pseudoscalar particles which is independent of the choice of the U(1)a invariant terms in Vq may be obtained 
by differentiating the following matrix equation representing the response of the potential to an infinitesimal axial 
transformation: 



dS 



dVo 
dcj) 



+ ((/., 5) -.(0', 5') = 1 



2Tr(0'||-5'||)-8c3^F,(A^,M') 



(13) 



To get general constraints on the pseudoscalar particle masses we differentiate this equation once with respect to each 
of the two matrix fields: (/>, (j)' and evaluate the equation in the ground state. Thus we also need the "minimum" 
condition, 



as 



SB 



as 



('^ 
\ as' 



avf 



SB 



OS' 



(14) 



III. MASSES AND MIXINGS 



As we previously discussed, the leading choice of terms corresponding to eight or fewer quark plus antiquark lines 
at each effective vertex reads: 

Vo^ - C2Tr(MMt) + c^Tr(MAf^MM^) 

+ d2 Tr(M'M't) + e'^ieabce'^'^MSM^M'/ + h.c.) 



C3 



, , detM , ,Tr(AfM't) 
.^^^"(d^) + (^-^^) Tr(Af-A^t) 



(15) 



All the terms except the last two have been chosen to also possess the U(1)a invariance. 
The minimum equations for this potential are: 



2a (-02 + 2040^ + 463/?) =2A, 



(16) 



avo 



2(d2/3 + 2e^a^) =0 



(17) 



Differentiating the potential in Eq. p5|) twice will yield four 2x2 mass matrices denoted as (Af^), (Afo), and 
[Xq) respectively for the pseudoscalar octets, the pseudoscalar singlets, the scalar octets and the scalar singlets. These 
may be brought to diagonal (hatted) form by the following 2x2 orthogonal transformations: 

Y,{R-^)ab{MI)bc{R.)cd = {MDad, J2(R^'^ab{Mo)bc{Ro)cd = {^Qad, 

B,C B,C 
^{La^)AB{Xj)Bc{La)cD {XI)aD, ^{L^^)ab{Xq)bc{Lo)cD = {Xo)aD, (18) 



B,C 



B.C 



Notice that the four mass matrices are identical to those given in section IV of [l| for the zero quark mass case. The 
numerical values of the entries will however differ because the relations among the coefficients are different due to the 
presence of 2 A rather than zero on the right hand side of Eq. ([Tl 
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FIG. 1: The predictions for the masses of the two SU(3) singlet scalars vs. m[7r(1300)]. The solid lines correspond to the 
massive pion case while the dashed lines correspond to the massless pion case previously considered. 

In the massive case there are 9 parameters (A, a, /?, C2, (i2, C4, ef, C3 and 71). These can be reduced to seven by 
use of the two minimum equations just given. We note that the parameters C3 and 71, associated with modeling the 
U(1)a anomaly, do not contribute to either the minimum equations or to the mass matrices of the particles which are 
not 0" singlets. Thus it is convenient to first determine the other five independent parameters. As the corresponding 
experimental inputs [3j we take the non-strange quantities: 



Evidently, a large experimental uncertainty appears in the mass of 7r(1300); we shall initially take the other masses 
as fixed at their central values and vary this mass in the indicated range. Essentially TO[7r(1300)] is being treated as 
an arbitrary parameter of our model. As shown in Eq. (|Aip in Appendix A, it is straightforward to determine the 
five independent parameters in terms of these five inputs. This determination is a generalization of the one in the 
previous zero mass pion case in which four parameters were determined from four inputs. 

The effects of adding a non zero quark mass term on the masses of the two predicted scalar singlets are displayed 
in Fig[TJ It is clear that the small mass term has a negligible effect on the mass of the heavier scalar singlet. On the 
other hand, there is a larger effect on the mass of the lighter scalar singlet. Still this singlet is exceptionally light so 
there is no qualitative difference in the result. 

In Fig. [2] we display a comparison of the four quark percentages of the tt meson, the lighter uq meson and the 
lighter scalar singlet with the corresponding values in the model with zero quark masses. (These are, of course, equal 
to the two quark percentages of the heavier particles with the same quantum numbers). It is clear that there is not 
much change compared to the zero quark mass case. 



It remains to discuss the four quark percentages of the two SU(3) singlet pseudoscalars. The lightest is the 7;(958) 
while candidates for the heavier one include ?7(1295), 77(1405), 77(1475) and 77(1760). As in the zero mass case, the 
first two candidates are ruled out because they do not lead to positive eigenvalues of the prediagonal squared mass 
matrix (Mq ). For the other two scenarios we may find the numerical values of the remaining parameters C3 and 71 
by using Eqs. (|A2|1 . (jA3p . (IA4[) and (|A5p given in the Appendix. The four quark contents for the 77(958) are being 
compared between the massive and massless quark cases in Fig. [3] Note that there are two solutions for each scenario 
corresponding to Eq. (|A4|) being of quadratic type. The lower four quark percentage curves seem the most plausible. 
Again there seems to be little difference between the zero and non-zero quark mass cases. This is understandable by 
comparing the values of the Lagrangian parameters found in Appendix A with those found in Appendix B of Jj] . 



77i(0''^octet) 
r7i(0+octet') 
r7^(0~octet') 

7n(0~octet) 



m[ao(980)] = 984.7 ± 1.2 MeV 
m[ao(1450)] = 1474 ± 19 MeV 
m[7r(1300)] = 1300 ± 100 MeV 
= 137 MeV 



131 MeV 



(19) 
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FIG. 2; Plot of the four quark percentages of various particles in the model as functions of the undetermined input parameter, 
m[7r(1300)]. Starting from the bottom and going up, the curves respectively show the four quark percentages of the pion, the 
lighter 0"'' singlet, and the ao(980). Solid lines apply for the massive pion case and dashed lines for the massless pion case. 
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FIG. 3; Plot of the four quark percentages of the r;(958) as functions of the undetermined input parameter, m[7r(1300)] for two 
scenarios. The top and bottom curves correspond to choosing the 7?(1760) as the heavier 0~ SU(3) singlet while the middle 
two curves correspond to choosing the 7^(1475) as the heavier 0~ SU(3) singlet. Note that for each scenario, the two curves are 
associated with different solutions of the quadratic equation (|A4[) for 71. The solid lines apply to the massive pion case while 
the dashed lines apply to the previous massless pion case. 



IV. PION SCATTERING: APPROXIMATE ANALYTIC TREATMENT FOR GENERAL Vo 



We have seen that there are not very big changes when using an SU(3) symmetric quark mass term of proper 
strength to give the experimental mass value to the pion. However, for the discussion of the pion pion scattering 
amplitude near threshold, which is one of the most important applications of the chiral approach, the correct value 
of the pion mass is important. In the zero quark mass case we showed that for any choice of terms in the potential 
Vq, the current algebra expression for the threshold pion pion amplitude held exactly. This was understandable since 
the algebra of chiral currents held by construction and furthermore the pion completely saturated the axial current. 
In the present case, the pion does not, as we shall see, completely saturate the axial current. Thus it is interesting to 
study this case in more detail. 

Note that the transformation between the diagonal fields (7r+ and 7r'+) and the original pion fields is given as: 



cos — sm f 
sin cos 6; 



(20) 



The value of the mixing angle, 6-,^ was previously written for an arbitrary potential Vq and with the symmetry breaker, 
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Eq.gl) as: 



wherein, 



tan(26'^) 



2A 

a 



d<t>'ldct>'l 



a 



(21) 



(22) 



The specific values of a;^, and will depend on the particular potential. In the case of no symmetry breaking, 
there was a big simplification and tan^^ was given by just -(3 /a. Furthermore, x^r and j/tt would be respectively the 
squared pion mass and the squared 7r(1300) mass in the absence of mixing. Clearly, the ratio x-^jy-r: is a very small 
number in the general case. We now make use of this fact to solve for the first correction to the mixing angle obtained 
from Eq. ([2T]l : 



sm( 



cos 



(a2+/?2) 



(a2 + f32y 



(23) 



In these equations the first terms correspond to the massless pion case and the second terms to the leading Xtt/j/tt 
corrections when the pion mass is turned on. 

An important application of Eq. (|23p is to the axial vector current with the pion's quantum numbers. After taking 
account of the mixing in Ea. ipPl) . the right hand side of the first Eg. pTI)) may be rewritten as: 



(24) 



where, 



= 2 a cos Ot, — 2 (3 sin 0.^ , 
F^/ = 2 a sin 6*^ + 2/3 cos 9.^ . 



(25) 



In the zero pion mass case, ^V' is seen to vanish. For the non zero pion mass case, we find, using Eq. (|23p . 



= 2V^ 



o\ — 



(a2 + p 



2^3/2 



o 



(26) 



It is seen that F-^ does not change much while Ftt' picks up a non zero value. Thus the tt' does not decouple from 
the axial vector current in the massive pion case. This means that PCAC does not strictly hold in the massive case 
and hence there is no reason to expect that the current algebra threshold theorem should be exactly correct in the 
present model. 

As a check of the accuracy of the approximate formula (j26p for F^/ we made an exact numerical calculation using 
the specific potential in Eq. p3|) and found F^' — —6.937 x 10^** GeV while, for the same parameters, the approximate 
formula gave F^' = -6.864 x 10"^ GeV. 

Now let us discuss the pion pion scattering in the threshold region. Our initial goal will be to see what results may 
be obtained for a general choice of chiral invariant^otential, Vq in Eq. ^ . The general pattern of this discussion and 
the notation is given in sections V, VI and Vll of [1] for the massless pion case. We start with the conventional tt — tt 
scattering amplitude at tree level: 



{XI) 



+ 



DD 



{XI 



(27) 



DD 



In this equation, g denotes the coefficient of the four point contact interaction among the physical (mass diagonal) 
pions. Furthermore g^D denotes the three point coupling constants connecting the physical pions to the two physical 
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SU(3) singlet scalar mesons. Similarly ggo stands for the coupling constants connecting the physical pions to the two 
physical scalar mesons which transform as the eighth component of an SU(3) octet. The usual Mandelstam variables, 
s,t, u are being employed. It is straightforward to numerically calculate the coupling constants just mentioned and 
the amplitude if the form of Vq is specified [e.g. Eq. (|15p ]: this will be discussed in the next section. To proceed with 
the general case we note that, for example, the coupling constant goo may be written as: 



Xdn+dn d{Sop)D/ 



di(l)l)Ad{(t>l)BdiSo)c 



(28) 



As discussed in fll| , there is a relation between the three point coupling constants on the right hand side and two point 
elements of the squared mass matrices. Such a relation follows from differentiation of Eq. (fT3|) together with the use 
of Eq. p4p . In fact it is similar in form to the relation obtained for the zero mass pion case: 



9oD = {Rt,)ai{Lo)hd [{XDah - [MDah] 



(29) 



wherein we have now adopted the convention of summing over repeated indices. The elements of the pion transfor- 
mation matrix, {R-„)ai are the angles sin0^ and cos^tt given in Eq. (|23p for the present non zero pion mass case. In 
the zero pion mass case, sin 0^ and cos 9-^ may be rewritten in terms of just a and /3 as we see by setting the second 
and higher order terms on the right hand sides of Eq. (|^5|) to zero. That is the form in which Eq. ([M)) exactly holds 
also for the massive pion case. Since Eq. (|29p and its analog for the four point vertices play an important role in the 
proof of the current algebra theorem, we can only prove the theorem in the massive case when we make the (not too 
bad) approximation that jy.^ is zero. 

With this approximation we can show, by generalizing the treatment given in that the usual "current algebra" 
formula holds for the massive pion case. In the previous treatment, the second term on the right hand side of Eq. ((29l) 
made no contribution. Now we must take this term's contribution into account. To get the same delicate cancellation 
due to chiral symmetry we must expand the amplitude in powers of s — instead of simply powers of s. Explicitly, 



A(s,t,u) = -- 



DD 



{Mi 



111 



DD 



(Mi) 




.981? 



90D 



{XDdd 



(M2) 



DD 



(30) 



Note that {X'^)dd is a single number indexed by D. There is a huge simplification of the coefficients of the (s — m^) 



term: 



5oD 



{XI 



dd 



^^(^7r^)lG(io)G£;(io ^)ea{XI)ah{Lo)hD 



(^TT ) Al ( Af^ ) AAf (^TT ) Af P (-Rtt ) -P-ff (^o) 



(M^)ii 
1 



(M)dd 

1 



-^^{RTj)Al{Ml)AM{RTi)Mp{RT^^)pH{Lii)HD 



{XDoD (M2)n 
1 



-2{RTr)ci[{XQ)cK — {M^)ck\{Lq)kd 
T2-(io '^)dk{XI)kc{Lq)cf{Lq ^)fj{R7t)ji 

11 

{Lq ^)dk{Xq)kc{Lo)cf{Lq '^)fj[Rtt)ji + 



{X^,)dd (M2) 
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•{R-iT ^)lc{M^)cN{Rpi)NF{RTr^)FK{Lo)Kl 



■^^iRT,^)lG{Lo)GE 



iXQ)DD 



{xi)dd {m)i 



-2 (-^0 )ej{R7,).J1 - 



(31) 
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Similarly, 



91d 



2 



(32) 



Next we must consider the contribution of the terms independent of (s — m^). The four point coupling constant g 
is approximately related to the matrix elements of the squared mass matrices (again with the proviso that the pion 
transformation matrices be similarly approximated) as: 



F2 



i(i?-i)iz, [{X^)dj - {M^)dj] {R.)ji + ^(i?;')ii? [{xS)dj - iM^)Dj] {R.).n 



(33) 



(The analogous calculation for the zero mass pion case is given in section VII of [1].) Using calculations similar to 
Ea. ([3T|) we also get: 



91d 



(Xi)nD-{M^)i 



(34) 



91d 



{XI) 



DD 



(M^)i 



llD 



Jl 



(35) 



Putting the last three equations into Eg. ([50]) we see that the sum of the terms independent of (s — m^) vanishes in 
the given approximation. The usual formula, 



Ais,t,u) = —{s-ml) 



(36) 



is thus obtained as an approximation. 

It should be remarked that Eq. p6p corresponds to keeping only terms up to linear order in the expansion for 
A{s,t,u). That means there are corrections, even at threshold, due to the masses of the scalars not being infinite. 
To see this and to summarize in a simple way, the preceding steps let us expand to one higher order, adopting a 
condensed notation in which rrii stands for the mass of any of the four scalars while gi stands for the corresponding 
trilinear coupling constant of that scalar with two pions. The first four terms in the expansion of A{s, t, u), as exactly 
obtained from Ea. ((27|) . are 



9 
2 



9t 



{s - ml) 



1 



+ (s - m^) ^ 



+ (- 



(37) 



The exact first equation contains, for each rrii, a geometrical expansion in the quantity (s — ml) /{ml — rnj). Thus 
the radius of convergence in s for this expression is the squared mass of the lightest scalar singlet. To apply this 
expression in the resonance region we must, of course, unitarize the formula in some way. Here we will be content 
to look at the threshold region. In going from the first to the second equation we used the facts established above 
that (in the approximation where F^' =0): 1) the sum of the first two terms of the first equation vanishes and 2) 
the third term of the first equation simplifies to becomes the first, current algebra, term of the second equation. The 
third term of the second equation represents the leading correction to the usual current algebra formula. It depends 
on the masses of the scalar mesons and would vanish in a hypothetical limit (often used) in which the scalar meson 
masses are taken to infinity. Note that every term in the approximate amplitude vanishes for s = m^, an unphysical 
point called the Adler zero [1]. Our derivation shows that the Adler zero follows from the generating equation (|13p . 
which in turn expresses the chiral invariance of the potential, Vq and from the saturation of the axial vector current by 
the pion field (so-called partial conservation of the axial current). The second equation in (j37p is an approximation, 
though a numerically good one, because the saturation of the axial current has been seen to be not strictly accurate 
in the present model. 

The situation in the case of zero pion mass [1,] is slightly different. There the amplitude is proportional to s so the 
Adler zero occurs at s = 0, which is also the threshold. Thus the current algebra amplitude as well as the corrections 
due to non-infinite mass scalar mesons vanish at threshold in the zero pion mass case. 

In the above we found that the current algebra theorem for a general potential does not seem to be exactly correct. 
This small deviation and in addition the more important effect of the scalar mesons will next be calculated exactly, 
by numerical means, for the scattering amplitude using the leading choice of Vq discussed in section III. 
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V. PION SCATTERING: EXACT NUMERICAL TREATMENT 



In the exact numerical treatment we do not need to make use of the relations between two point and three point 
functions and between three point and four point functions since we adopt the specific potential Vq given in Eq. 
P^ . The needed quantities for calculating the scattering amplitude are displayed in Eq. (|77)): the four physical 
scalar singlet masses, the four three-point coupling constants conecting these scalars to two pions and the four-point 
pion physical coupling constant, g. These are obtained by, in turn, differentiating the potential twice with respect 
to two scalar fields ( gf gg ) , two pseudoscalar fields as well as one scalar field ( gsd<i>d<p ) pseudoscalar fields 

( QjjOd>dibdd) ) • Furthermore we must use equations like Eq. (^5)) to relate the "bare" amplitudes obtained by such 



differentiations to the physical ones (ie, in mass diagonal bases for the fields). The matrices transforming the fields 
to mass diagonal bases are defined in Eq. (jlSp and are obtained by diagonalizing the relevant squared mass matrices. 

For our purpose we define the current algebra result in terms of the expansion of the tree level amplitude A(s, t, u) 
in powers of (s — mj), as displayed in Eq. ([57)1 . Specifically, 



(38) 



where, 



Co 



E 



9i 



(mf - m2 )2 ■ 



(39) 



The current algebra result requires Co to vanish and Ci = S/F^. 

Plots of Co and Ci as functions of the model parameter m[7r(1300)] are shown in Figs. 2] and O respectively. Even 
though Co is small it is clearly non vanishing. Also Ci deviates by a few percent from the current algebra prediction. 
To estimate the numerical accuracy of this calculation it was repeated for the case of zero pion mass. There it was 
found that Cq = 0(10^*) whereas it should be exactly zero. Thus the accuracy of the calculation method is several 
orders of magnitude more sensitive than the indicated effect. In this model, the Adler zero is shifted (by about 
-Co/Ci) very slightly to the left of . 




1.22 1.23 
111 11(1300)1 (GeV) 



FIG. 4: Constant term, Co in expansion of invariant amplitude as a function of the model parameter m[7r(1300)]. The error 
bars reflect the uncertainty in the mass of ao(1450). Current algebra would predict Co=0. 

The small deviations from the current algebra result just discussed seem to be beyond present experimental accuracy. 
On the other hand, the "beyond" current algebra contributions to A{s,t,u) due to the higher than linear terms in 
the expansion shown in Eq. (|37|) seem to be highly relevant for comparison with present day experiments 0, @, S S] ■ 
These corrections would vanish in a limit where the scalar masses all go to infinity, which essentially corresponds to 
the use of a non-linear rather than the present linear type of sigma model. 
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1.22 1.23 
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FIG. 5: Coefficient of linear term Ci (in units of GeV ^) in expansion of invariant amplitude compared with the current algebra 
result (dashed line) as a function of the model parameter m[7r(1300)]. The error bars reflect the uncertainty in the mass of 
ao(1450). 



It is usual to discuss the amplitudes near threshold in terms of their partial wave scattering lengths. The J = 
scattering lengths are of course especially affected by the presence of light scalar mesons. Using the compact notation 
in Eq. ([57|) . the (dimensionless) partial wave scattering lengths may be calculated to be: 



rriTra, 



1 
1 

32^ 



5ff , 2 



^ mf 



3 



(40) 



These formulas are expressed in terms of the physical masses and coupling constants which are being computed 
exactly by numerical means. Note that the isospin label, / and the angular momentum label, J appear as Oj. For 
comparison, we may give the usual current algebra results : 



7mi 
167rF2 

-2ml 
167rF2 ■ 



(41) 



The results of our numerical calculation are shown in Fig. [6l 

It is seen that the numerical calculation for the scattering length in the non resonant 1 — 2 channel gives about 
the same value, -0.04 as the current algebra result. In the resonant 1 = s-wave channel the current algebra result of 
0.15 is smaller than the result of the exact calculation for the range shown of the model parameter TO[7r(1300)]. The 
exact calculation result in this resonant channel varies strongly with m[7r(1300)] in contrast to the case in the non 
resonant channel. To understand what this means we should ask what is the significance of varying m[7r(1300)]. Such 
a variation might be associated with variations in the masses of the four iso-singlet scalars. But the set up of the model 
as shown in Eq. (fT9| fixes the physical masses of the two octet isosinglets (at 985 MeV and 1474 MeV) . Furthermore 
we see in Fig. [T]that varying r7i[7r(1300)] leaves the mass of the heavier SU(3) singlet scalar essentially unchanged 
(at about 1500 MeV) while it changes the mass of the lightest SU(3) singlet scalar. Thus it seems unavoidable to 
interpret the variation of r7i[7r(1300)] as being associated with the variation of the mass of the lightest scalar. This is 
confirmed by noticing that the largest change in aj] occurs in the region of TO[7r(1300)] where the mass of the lightest 
scalar is changing most rapidly. 

The correction to the current algebra result for Gq due to the finite masses of light scalars was already discussed 
and noted to be positive in [oj, some years ago. The contribution of a light scalar meson to the scattering length was 
recently calculated in [Toj . 
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1.22 1.23 
m lll(1300)J (GeV) 



FIG. 6; Top curve: I — J = scattering length, m^aSj vs. m[7r(1300)]. Bottom curve: I — 2, J = scattering length, m.n-ao 
vs. m[7r(1300)]. The error bars reflect the uncertainty of m[ao(1450)]. 



DIRAC Collaboration 



It is very interesting to examine the recent experimental data on the s-wave scattering lengths Oq and Qq] these 
include the following. 

NA48/2 collaboration Q: 

m^+(a[!-a^) = 0.264 ±0.015 (42) 
m^+a[] = 0.256 ±0.011 (43) 

E865 Collaboration [1]: 

m^+aj] = 0.216 ±0.015 (44) 

m,+ ag = 0.264;;^;°38 (45) 
A general discussion of these experiments is given in Q . 

Comparison of experiment with theory shows that the larger values of Oq predicted by the numerical calculation 
when m[7r(1300)] is greater than about 1215 MeV, give good agreement. (This corresponds to the lightest scalar singlet 
lighter than about 460 MeV). In contrast, the current algebra prediction for Oq is clearly too low. The nonresonant 
channel with (/ = 2, J = 0) is not so well determined from experiment but seems to be consistent with the common 
prediction of current algebra or the numerical calculation. 

The indicated value of the lightest scalar mass in our model is consistent with recent results [ll| obtained by using 
Roy dispersion relation sum rules. The typical values obtained for the mass and width of the lightest scalar are M 
= 441 MeV and T = 544 MeV. These are also similar to what is obtained jj], M = 457 MeV and T = 632 MeV, 
by using a K-matrix unitarized three flavor linear sigma model (with just one chiral nonet). A unitarized two flavor 
linear sigma model was earlier given in (l3| . The s wave pion pion interaction has recently [l^ been discussed using 



the Adler sum rule. Actually, a long time ago the Adler sum rule was used [15[ to suggest a light scalar with a similar 
mass to the above. 

The encouraging result for Qq with inclusion of scalar meson corrections corresponds to a tree level treatment of this 
linear sigma model. One may justifiably wonder whether the agreement would be spoiled by inclusion of loops, i.e. by 
a unitarization of the model. While this is a more complicated matter we may note that the simplest unitarization, 
the K matrix approach, does not change the result at all. In this approach the corrected scattering length would be, 

m^ajj (corrected) = ■ —=====, (46) 

1 — tniT^aQy'l — Am^/s 

where the kinematical square root is to be evaluated at threshold, s = 4m^. Of course this is just an indication 
rather than a proof that the effect of unitarization is expected to be small for the scattering lengths. 
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FIG. 7: Individual contributions in Eq. (|27p to A{s,t,u) at threshold. 
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FIG. 8: Individual contributions to the second of Ea. (|37p to A{s,t,u) at threshold. 

VI. DISCUSSION AND CONCLUSIONS 

We have seen that the intuitive explanation for the existence of a very Ught scalar meson as well as light scalar 
mesons vifith large "four quark" content given in [l[ (which used the simplifying assumption of zero pion mass) still 
is good when the physical pion mass is used. The main physical input is (almost) spontaneously broken chiral 
symmetry. Our model treats the chiral symmetry in its linear realization. This is normally considered inconvenient 
for the treatment of the scattering problem since it is known that there are very large cancellations between the four 
pion contact and scalar meson pole terms due to the symmetry. However it has the nice feature that it retains the 
light scalar mesons which give crucial corrections to the current algebra results. The inconvenience due to the large 
cancellations can be removed by using the Taylor expansion in Eq. ([37]). The third term in the second line of that 
equation shows that the beyond current algebra corrections can be neatly displayed in a physical way; they are seen 
to be of order [ni.,^/niiY compared to the current algebra term. Here is a scalar meson mass. Still higher order 
corrections are suppressed by additional orders of {nriT^/mi)'^. 

It is very interesting to illustrate how the use of the expanded form of the amplitude given in Eq. ([57]) leads to a 
nice picture. In Fig. [7] the individual contributions to the starting form of the amplitude, Eq. (|27p . associated with 
each scalar exchange as well as with the four point contact term arc displayed. Evidently the contact term is the 
largest in magnitude and cancels off most of the other contributions. In addition the contribution from the lightest 
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scalar exchange, using this starting equation, has a very small magnitude, although its effect is expected to be 
the largest. Evidently, the intricate cancellations completely distort the underlying physics when expressed in this 
form. On the other hand, the use of the second of Ea. ([57|) clears things up, as may be seen from Fig. [S] There the 
correction to the current algebra result is seen to be completely dominated by the lightest scalar. 

Clearly the s-wave pion pion scattering is rather nicely described by a linear sigma model. At the present stage it 
does not seem to matter which one is chosen for this aspect. The virtue of the present toy model is that it accomodates 
both 2 quark and 4 quark scalar nonets in a consistent way and may help in understanding the relation to QCD as 
earlier discussed [l[ . Another interesting consequence of this linear formulation is the presence of heavy, largely four 
quark pseudoscalars, of which the 7r(1300) is a possible candidate. 

An amusing feature of the present model with a massive pion is that, due to a small deviation from exact PCAC, 
the current algebra result for near threshold pion pion scattering no longer holds exactly. Associated with this is 
the feature that the amplitude no longer vanishes at the unphysical point (Adler zero), s = m^. The zero is shifted 
somewhat. The measure of this effect is the deviation of the coefficient, F^i from zero. We have seen that this is 
quite small since F^^i is of the order x^/y^ ~ m^/m^, ~ 0.01, as may be seen from Eq. (pS)) . On the other hand the 
similar effect should be more noticeable for the K — K' system and ior K — K scattering. In [l^ it is shown that this 
system is formally analogous to the tt — tt' system with the substitutions x-^; — > xk and y,r Uk- But in this case 
Fk' is expected to be of the order rn\/m\, « 0.1, ten times greater than for the pion case. A full calculation would 
require setting up the model with the inclusion of SU(3) symmetry breaking quark mass terms. This will be studied 
elsewhere. 

A possible general question about the present model is that it introduces both states made of a quark and an 
antiquark as well as states with two quarks and two antiquarks. According to the usual 't Hooft large Nc extrapolation 
[l7l | of QCD the "four quark" states are expected to be suppressed. However, it was recently pointed out [Ts] that the 
alternative, mathematically allowed, Corrigan Ramond [l9| extrapolation does not suppress the multiquark states. 
This kind of extrapolation may be relevant for understanding the physics of the light scalar mesons. 
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APPENDIX A: PARAMETER DETERMINATION 

Given the inputs: the pion decay constant, F^^] the mass of the pion, m^; the mass of the ao(980), ma] the mass of 
the ao(1450), nia'] the mass of the 7r(1300), m^r', the independent model parameters which don't involve the U{1)a 
violating terms can be successively determined (in the order given) by the equations: 
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Once the above parameters are determined, the parameters 71 and C3 of the U{X)a violating sector are obtained 
in terms of the mass of the 77(958), m^i and the mass of a suitable heavier 0~ isosinglet, m^2 using the following 
procedure. The 2x2 prediagonal mass-squared matrix of the two SU(3) singlet pseudoscalars is written in the form: 
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(A2) 



where Kij is a real symmetric matrix involving the coefficients of the terms in Vq which are U(1)a invariant . With 



the choice of invariant terms in Eq. (|T5|) we have 

Kn = 
K12 = 
K22 = 



-8 eg a 



2 do. 



(A3) 



Then, 71 is found as a solution of the quadratic equation: 
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In addition, 
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Next we give the numerical values of the parameters for the central values of all the listed input masses except for 
m[7r(1300)] which instead will take the typical value allowed by both the data and by the model, 1215 MeV. Table |I] 
shows the results for the parameters which are not associated with the U(1)a violating part of the Lagrangian. 



C2(GeV2) 


8.79 


d2(GeV2) 


6.30 xlO"^ 


eg(GeV) 


-2.13 




42.4 


a(GeV) 


6.06 xlO-2 


/3(GeV) 


2.49 xlO"^ 


A(GeV^) 


6.66 xlO"'' 



TABLE I: Calculated Lagrangian parameters:c2, d2, 63, cf and vacuum values: a, (3. 



Table HIl shows the calculated Lagrangian parameters associated with the U(1)a violating terms. Two "scenarios" 
associated with different identifications of the heavy rj which is the partner of the 7y(958) are shown (I assumes ?7(1475) 
to be chosen while II assumes ri{1760) to be chosen.) For each scenario, the two solutions (labeled 1 and 2) are shown. 

Using these parameters we next list the mixing matrices for, respectively, the two 0~ octet states, the two O"*" octet 
states and the two 0"*" singlet states: 
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11 


12 


III 


112 


C3(GeV^) 


-2.39 X 10"* 


-2.38 X 10"* 


-3.42 X 10~* 


-3.37 X 10"* 


71 


5.33 xlO"^ 


2.52 X 10"^ 


8.68 xlO"^ 


-8.65 X 10-2 



TABLE 11: Calculated parameters; cs and 71. 



0.923 0.385 
-0.385 0.923 



-0.486 0.874 
0.874 0.486 



0.706 0.708 
0.708 -0.706 



Similarly, the mixing matrices for the two solutions for scenario I of the singlet states are: 



I 1 : (Ro') 



-0.671 0.742 
0.742 0.671 



, /2:(i?oi) 



0.853 -0.522 
0.522 0.853 



Finally, the mixing matrices for the two solutions for scenario II of the singlet states are: 



(A6) 



(A7) 



-0.411 0.912 
0.912 0.411 



II 2 : {Ro') - 



0.972 -0.235 
0.235 0.972 



(A8) 
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